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We calculate the effective fluctuation induced force between spherical or disk-like colloids trapped 
at a flat, fluid interface mediated by thermally excited capillary waves. This Casimir type force is 
determined by the partition function of the system which in turn is calculated in a functional integral 
approach, where the restrictions on the capillary waves imposed by the colloids are incorporated 
by auxiliary fields. In the long-range regime the fluctuation induced force is shown to depend 
ly-^ sensitively on the boundary conditions imposed at the three-phase contact line between the colloids 

^vq and the two fluid phases. The splitting of the fluctuating capillary wave field into a mean-field and a 

fluctuation part leads to competing repulsive and attractive contributions, respectively, which give 
rise to cancellations of the leading terms. In a second approach based on multipole expansion of the 
Casimir interaction, these cancellations can be understood from the vanishing of certain multipole 
moments enforced by the boundary conditions. We also discuss the connection of the different types 
of boundary conditions to certain external fields acting on the colloids which appear to be realizable 
■ by experimental techniques such as the laser tweezer method. 
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PACS numbers: 82.70.Dd, 68.03.Kn 



I. INTRODUCTION 



The structure formation of nanoscopic colloidal particles adsorbed at fluid interfaces has attracted an in- 
creasing interest in recent years because of the various applications such systems exhibit, as in the design of 
, nanoscale devices 1] (e.g. for optical applications), and also because of the physical insights offered, e.g. for 
£" — ■ the understanding of protein aggregation on cell membranes. The practical importance arises from the fact 
lf~) \ that nano- and microcolloids are very effectively trapped by fluid interfaces [2(. The high stability of partially 
wetting colloids (with their size ranging from nm to /jm) at interfaces enables the formation of two-dimensional 
ordered structures and also of rather complex mesoscale patterns (c/. Refs. |E 0, IE IE IE IE EH E])- in 
spite of the numerous experimental and theoretical efforts in the last decade, the nature of the effective forces 
between the colloids governing the arrangement of the colloids trapped at the interface is not fully understood 
yet. This holds in particular for the mesoscale pattern formation which points to sizeable and quite long-ranged 
attractions between the colloids. For the colloidal interactions at interfaces three different regimes may be 
distinguished. For colloid sizes of about 0.5 — 5 /im, capillary forces are important. Though gravity is known to 
I ■ be unimportant in this regime, electrostatic forces (either caused only by the charge distributions on the colloid 
surface and the fluid phases, or additionally impose d by an external field) may lead to interfacial deformations 
and, hence, to considerable capillary interactions 0,E3>E1- I* 1 this regime, also strongeffects of colloid surface 
inhomogeneities leading to an undulated three-phae contact line appear to be relevant |j| ■ In the opposite limit 
of colloid sizes of a few nanometers, effects of interparticle correlations within the interface become important 
and determine the structural properties of the system. In this regime the colloidal interactions should be treated 
by truly microscopic fluid theories like density functional or inhomogeneous integral equation approaches (for 
a general scheme, see Ref. EH> an d f° r results on correlations within a free interface, see Refs. EEE3) In 
between these two regimes, the microscopic one which must be tackled with the full power of classical statistical 
mechanics and the macroscopic one where thermal fluctuations appear to be unimportant, a coarse-grained pic- 
ture of the fluctuating fluid interface should be applied. Within such a picture, the properties of fluid interfaces 
are very well described by an effective capillary wave Hamiltonian which governs both the equilibrium interface 
configuration and the thermal fluctuations (capillary waves) around this equilibrium (or mean-field) position. 
As postulated by the Goldstone theorem the capillary waves are long-range correlated. The interface breaks the 
continuous translational symmetry of the system, and in the limit of vanishing external fields - like gravity - it 
has to be accompanied by easily excitable long wavelength (Goldstone) modes - precisely the capillary waves 
EH- The fluctuation spectrum of the capillary waves will be modified by colloids trapped at the interface and 
therefore leads to fluctuation-induced forces between them. In the case of anisotropic colloids (rods) these forces 
have been evaluated in Ref. E3 an d shown to lead to an orientational dependence. Furthermore we note that 
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there is numerous work on the force between inclusions on membranes where the membrane shape fluctuations 
take the role of capillary waves , see e.g. Refs. In this paper, we shall address the fluctuation-induced 

effects of capillary waves on rotationally symmetric colloids, i.e., spherical or disk- like ones, and pay special 
attention to the effects of the boundary conditions at the three-phase contact line (where the fluid interface 
meets the colloid surface) on the large-distance behavior of the fluctuation-induced force. In previous work [2l| 
we have established the independence of the short-distance behavior on the boundary conditions. For colloids 
at contact, the fluctuation-induced force is attractive and diverging, albeit the div erge nce is somewhat slower 
than for the ubiquitous van der Waals (vdW) forces. However, as discussed in Ref. |2l[. the fluctation-induced 
force dominates for specially prepared system (such as index-matched ones) and hence it leads to colloidal 
aggregation at the interface. 

Put in broader context, the prototype of such a fluctuation induced force as considered here is the original 
Casimir effect of the zero-temperature, long-ranged quantum fluctuations of the electromagnetic field confined 
between two conducting plates [22j that was predicted by Hendrik Casimir in f 948. In f 978 this was extended by 
Fisher and de Gennes [22j to critical systems like bulk fluids near the critical point or in binary mixtures near 
their critical demixing point where instead of long-ranged quantum fluctuations one has long-ranged critical 
fluctuations |24| . An even more common example for long-range correlated fluctuations in thermal systems 
are the already mentioned Goldstone modes accompanying broken continuous symmetries. These modes are 
ubiqitous in condensed matter physics with famous examples being director fluctuations in liquid crystals 25], 
phase fluctuations in superfluids |26j or capillary waves at fluid interfaces mentioned above |27| . 

Both, the great experimental progress (see e.g. Ref. J28|) and the fundamental interest in the Casimir effect 
have triggered numerous studies in this field in the last twenty years. In spite of this effort, still the only 
exactly solvable case is Casimir's original one of two parallel infinite planes. An almost closed solution has been 
obtained for the plate-cylinder case very recently |29j. In more complex geometries, analytical results usually 
are only available in asymptotic ranges, and the strong dependence of the Casimir force on the geometrical 
arrangement of the system has been shown in the case of the fluctuating medium being a free scalar field for 
various configurations like, e.g., immersed metallic spheres [3(| or cylinders [SjJ and plates. A standard, though 
uncontrolled, approximation for the Casimir interaction between one curved body and a plate or two curved 
bodies at short distances is the well-known proximity force approximation which presumes pairwise additivity 
and which is in the same spirit as the Derjaguin approximation |3l| well known from calculations of the van der 
Waals force between macroscopic bodies. The limitations of this approach are, e.g., investigated perturbatively 
and numerically in Refs. [32Ll33| . 

In a study on the pairwise addititvity of the Casimir force |34| Golestanian investigated the force between 
spheres in a fluctuating medium described by a free scalar field in arbitrary dimension D. Using a terminology 
borrowed from electrostatics, the scalar "field" corresponds to the potential and if the spheres are presumed to 
be metallic, one can distinguish two types of boundary condition: (i) the spheres are grounded, corresponding 
to zero potential at their surface (Dirichlet boundary conditions) and (ii) the spheres are isolated with constant 
charge, in which case one averages over the surface potential. The asymptotic Casimir force between two 
spheres at distance d strongly depends on the type of boundary condition; in case (i) the author finds -Fground ~ 
l/d 2 (- D ~ 2 ) +1 and for case (ii) F mo ~ l/d 2D+1 . In relation to this study, our present work is only concerned with 
the case D = 2. The boundary conditions introduced above correspond to the physically realizable conditions 
of a fixed colloid with pinned three-phase contact line (grounded case) and a verticallly fluctuating colloid 
with pinned contact line (isolated case). Furthermore we will identify two more physically realizable boundary 
conditions (arbitrarily fluctuating colloid with pinned contact line and colloid with freely fluctuating contact 
line) which lead to different results for the Casimir force again. Also we will find that in D = 2 the general 
result for F groun d is modified since care is needed to perform the limit to a massless fluctuating field. We 
will elucidate the appearance of the leading terms in an asymptotic expansion of the Casimir force from two 
different computational schemes. In the first one, the force is calculated by splitting the fluctuating field into a 
mean field and a fluctuation part as, e.g., it is well known from the functional integrals in quantum mechanics. 
Then the partition function consists of two parts which can be computed separately and lead to a repulsive 
and an attractive contribution to the Casimir interaction, respectively. The leading non-vanishing term in the 
asymptotic behaviour for large colloid separations then strongly depends on the type of boundary conditions 
under consideration. In an alternative approach this can be understood from an auxiliary multipole expansion 
of the interaction between the colloids similar as in Ref. [34| . One can relate each type of boundary conditions 
to a suppression of certain auxiliary multipolcs on the colloids and the leading term of the Casimir force is 
obtained from the first non-vanishing multipole-multipolc interaction. 

The paper is organized as follows. In Sec. [n] we will introduce the model system by an effective Hamiltonian. 
We will show how the partition function can be calculated via functional integrals in two distinct ways mentioned 
above and how the various boundary conditions can be implemented. Then we will compute the corresponding 
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functional integrals in Sees. Hill IIVI and Ivl in the asymptotic regimes of long and short colloid separations d and 
numerically in the full range of d, respectively, and compare analytical with numerical results in the discussion 
part, Sec. I VII In Sec. I VII we will also discuss possible experimental tests for this type of Casimir force through 
the implications of an additional external potential for the colloids, which may, e.g., be realized by optical 
tweezers. 

II. MODEL 

In this section we derive an effective Hamiltonian for the free energy changes associated with thermally excited 
height fluctuations of the interface between two fluid phases I and II at which two nano- or microscopic colloids 
are trapped. As described above, this configuration is very stable against thermal fluctuations for colloids with 
a partially wetting surface. The colloids are assumed to be either spherical with radius R or disklike with 
radius R and thickness H . In the absence of charges and for colloid sizes R < 1 /im the weight of the particles 
can be neglected. Thus the equilibrium configuration of minimal free energy is the flat interface and spherical 
colloids are positioned such that Young's law holds at the horizontal three-phase contact line which is a circle 
with radius r = R sinO. Young's angle, measured through phase II (assumed to be of higher density than 
phase I), is determined by cos6> = (71 — 711V7 where 7 is the surface tension of the interface between I /II, 
and 7i[n] is the surface tension between the colloid and phase I [II], respectively (c.f. fig. For disklikc 
colloids, the contact line is either the upper (6 < tt/2) or lower (9 > 7r/2) circular edge, so that its radius is 
given by ro = R. For both spheres and disks the cross-section of the colloids with the flat interface is given 
by the interior of a circle which we refer to as Si^ e f below, and, therefore both cases can be treated within 
the same model in the following. We take the flat interface S'mcn,rof as the reference configuration with respect 
to which free energy changes are measured, and choose it to be the plane z = with the two circular holes 
S^ref of radius ro in it: S'mcn.rcf = R 2 \ UiS^ref (see fig. EJl- Deviations from the planar reference interface 
z = are assumed to be small and without overhangs or bubbles. That allows using the Monge representation 
(x, y,z = u(x, y)) = (x, z = it(x)) as a parametrization of the actual interface positions. The free energy costs 
for thermal fluctuations around the flat reference interface is determined by the the change in interfacial energy 
of all interfaces (I /II, colloid/I and colloid/II): 

Wtot = 7A^mc„ + J1AA1 + tqAAq . (1) 

The first term in Eq. Q expresses the energy needed for creating the additional meniscus area associated with 
the height fluctuations and is given by 

7AA mcn = 7 / d 2 iVl+(Vti) 2 -7 / d 2 x 

JS men "'Smcn.rcf 

d 2 x ( Vw) 2 + 7A^ P roj ■ (2) 

In Eq. P)l. S meTL is the meniscus area projected onto the plane z = (where the reference interface is located) and 
•Smen.ref is the meniscus in the reference configuration mentioned above. AA pro j = f s ^ s drx is the change 
in projected meniscus area with respect to the reference configuration. In the second line we have applied a 
small gradient expansion which is valid for slopes Vu <C 1 and which provides the long wavelength description 
of the interface fluctuations we are interested in. Note that the first line of Eq. @ constitutes the drumhead 
model which is well-known in the rcnormalistion group analysis of interface problems, but is also used for the 
description of elastic surfaces (c.f. Ref. |l8|). 

As discussed in the Introduction, the center of colloid i may fluctuate vertically (measured by hi) around 
the reference position as well as the contact line itself may do. This leads to changes in the interfacial areas 
colloid/I[II] (AAipij / in Eq. (Q). In order to determine the corresponding energy costs we introduce the 
vertical position of the contact line at colloid i as a function of the polar angle defined on the reference 
contact line circles 9Si, r ef) by its Fourier expansion 

00 

h=u{dS^ ci )= ]T P im e im ^ (3) 

m— — 00 

and refer to the Fourier coefficients Pi m as boundary multipole moments below. Since /, is real, we have 
Pim — P*- m - Following Ref. 0, the free energy changes associated with AA pro j and AAjpj can be expanded 




FIG. 1: Side view on the reference configuration (here the colloids are assumed to be spheres). 



up to second order in hi and /j, and may therefore be collected in a boundary Hamiltonian TCi^ (see App. lA")l : 

H h ,i[fi, hi] = 7AA proj + 71AA1 + 711AA11 
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2(-Pjo — fH 



m>l 
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Putting Eqs. Q,J2J an d together, the total free energy change is the sum 



Wb,i + Hb- 



d 2 x (Vu) 2 + H h ,x + n h , 2 



(4) 



(5) 



of the capillary wave Hamiltonian Tt cw which describes the free energy differences associated with the additional 
interfacial area created by the height fluctuations, and the boundary Hamiltonians Hb,i combining all effects 
related to fluctuations of the three-phase contact line on the colloid surfaces. As is well-known, the Hamiltonian 
Hew is plagued with both a short-wavelength and a long-wavelength divergence which, however, can be treated 
by physical cutoffs. The natural short- wavelength cut-off is set by the molecular length-scale a of the fluid at 
which the capillary wave model ceases to remain valid. The long wavelength divergence is reminiscent to the 
fact that the capillary waves are Goldstone modes. Of course, in real systems the gravitational field provides 
a natural damping for capillary waves. Accounting also for the costs in gravitational energy associated with 
the interface height fluctuations, therefore, introduces a long wavelength cutoff and leads to an additional term 
("mass term") in the capillary wave Hamiltonian, 

,2 - 



d z x 



(Vu) 
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a! 



(0) 



Here the capillary length is given by A c = [7/dpn — Pi\ s)] 1 ^ 2 ; where p t is the mass density in phase i and g is 
the gravitational constant. Usually, in simple liquids, A c is in the range of millimeters and, therefore, is by far 
the longest length scale in the system. In fact, here it plays the role of a long wavelength cutoff of the capillary 
wave Hamiltonian 7i cw , and we will discuss our results in the limit A c 3> R and A c 3> d. However, as we will see 
below, care is required when taking the limit A c — ► 00 (corresponding to g — > 0), since logarithmic divergencies 
appear [35| . Another common way to introduce a long- wavelength cut-off is the finite size L of any real system. 
As discussed in Ref. the precise way of incorporating the long-wavelength cut-off is unimportant for the 
effects on the colloidal length scale. As an example, in both approaches the width of the interface related to the 
capillary wave is logarithmically divergent, (u(0) 2 ) ~ lnX c [L]/a. 

Via the integration domain of Ti. cw , the total Hamiltonian of the system, Eq. JSJ, implicitly depends on the 
geometric configuration. This leads to a free energy which depends on the mean distance d between the colloid 
centers and gives rise to an effective force F(d) = —(dJ-)/(dd) as function of the mean local distance between 
the colloid centers which is determined by the free energy J-(d) — — kBT\nZ(d). The partition function Z{d) is 
obtained by a functional integral over all possible interface configurations u and fi\ the boundary configurations 
are included by ^-function constraints, 



Z' 1 
^0 



Vu exp 



n 



XiGdSi, rB f 



/i(xj)]exp-l - 



'Hb ,j\Jh hj 

k B T 



(7) 



-> 00) = I and ensures a proper regularization of the functional 
u is coupled to the contact line height /j and therefore, the 



Here Zq is a normalization factor such that Z[d - 
integral. Via the <5-functions the interface field 
boundary Hamiltonians Hi,b have a crucial influence on the resulting effective interaction between the colloids 
as we will see below. In the next subsection we discuss possible situations for the boundary conditions at the 
three phase contact line and specify the corresponding integration measure Dfi. 



5 



y 

d 




FIG. 2: Reference configuration with two circles Si tls f representing the reference contact line on the colloid surfaces 
(here the colloids may be disks or spheres). 

A. Boundary conditions at the three-phase contact line 

We shall discuss two different realizations of the boundary conditions for the contact line, namely 

• Case (A): The contact lines and the vertical colloid positions fluctuate freely; this corresponds to the 
physical situation of smooth, spherical colloids. In this case, the integration measure is given by T>fi — 
J dhi Y[ m dPim and encompasses integration over all boundary multipolc moments. 

• Case (B): The contact lines (the circles dSi tie i in the reference configuration) are pinned to the colloid 
surface. This corresponds to disklike colloids or Janus spheres consisting of two different materials, or 
to colloids with a very rough surface. Within the pinning case, we furthermore distinguish the following 
three situations: 

(Bl) The colloid positions are frozen, thus there are no integrations over the boundary terms. 

(B2) The vertical positions of the colloids fluctuate freely, thus boundary monopoles must be included in 
the integration measure so that X>/, = dhi dPio 5{Pio — hi). 

(B3) The vertical position and the orientation of the colloids (tilts) fluctuate freely. Up to second order 
in the tilts this corresponds to the inclusion of boundary dipoles in the integration measure, thus 
Vfi = dhi dP l0 dP a dPi-x 6{P l0 -hi), 

The (5-function expresses the pinning condition. 

In the following subsections we introduce two distinct schemes to compute the partition function Z(d) in 
Eq. (J7J which allows us to discuss the results from different perspectives. The first one is based on a splitting of 
Z in a mean- field and a fluctuation part, whereas for the second one we extend the fluid interface artificially 
- to the interior of the reference circles Si fTe f. As we will see, in the first approach, the mean-field part leads 
to a repulsive and the fluctuation part to an attractive Casimir force. The interplay of these two contributions 
gives rise to an cancellation of the leading terms in the asymptotic range ro <C d which sensitively depends 
on the chosen boundary conditions. In the second approach, the asymptotic Casimir force can be expanded 
in multipole interactions of auxiliary fields defined on S^ref & n d the appearance of a different behaviour in 
the long-range regime for different boundary conditions can be understood straightforwardly as the result of 
vanishing multipolc moments of the auxiliary fields, enforced by the boundary conditions. 
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B. Mean- field and fluctuation part 



As the capillary wave Hamiltonian is Gaussian in the field u of the local interface position, a standard 
procedure for the evaluation of the functional integral Q is the decomposition into a mean-field and a fluctuation 
part, 

U = Umf + V . (8) 

The mean- field part solves the Euler-Lagrange equation of the capillary wave Hamiltonian 7Y CW , 

(-A + \- 2 ) u„a = (9) 

with the boundary condition u m { \ds ilc{ — /•»• Consequently the fluctuation part is pinned to zero at the 
contact line, v rof = 0. Then the partition sum Z = Z nuc Z m t separates into a product of a fluctuation part 
independent of the boundary conditions, and a mean field part depending on the boundary conditions (which 
may fluctuate theirselves, see the cases (A), (B2), and (B3)): 

-2fluc = Zq 1 I T>v Yl J\ <K w ( x i)) ex P 
i=ix«e3S'i, r rf 

Z m{ = fljvfiexpL^^J^ ^/ i (x i )(5 nUmf (x i ))|exp|-^0M} . (10) 

The first exponential in Z m { stems from applying Gauss' theorem to the energy associated with u m f. In this 
term d n u m { denotes the normal derivative of the mean-field solution towards the interior of the circle dSi^ c {. 

For intermediate asymptotic distances between the colloids (ro <C d <C A c ), Zfi uc and Z m { are analyzed 
separately in Sees. IIII Al and UlI Bl respectively. 

This procedure to split the functional integral into a mean field part and a fluctuation part with a zero 
boundary condition is equivalent to the standard calculation of the quantum mechanical path integral for 
a single particle in a harmonic potential. The mean field part corresponds to the classical action and the 
fluctuation part evaluated becomes the fluctuation determinant (see, e.g., Ref. [37j ). 




C. Alternative approach (Kardar's method) 

An alternative scheme to calculate the partition function of the two colloids and the fluctuating interface 
without splitting it into a mean-field and fluctuation part can be devised if the interface height field u(x,y) 
which enters the functional integral for Z is extended to the interior of the circles S^ref- Thus the measure 
of the functional integral for Z is extended by 2?zt(x)| x6 5 i rof and the integration domain in the capillary wave 
Hamiltonian is enlarged to encompass the whole R 2 . We note that physically the free energy of the system 
must not change since in the interior of Si lIe i the interface is pinned to the colloid surface. A method similar 
to this ansatz was introduced in Ref. [19| by Kardar et al. which investigates effective forces between rods on 
fluctuating membranes and films, and therefore we refer to it as Kardar's method in the following. 

On the colloid surfaces, the interface height field is given by the three phase contact line, u(dSi tre {) = fa. We 
extend u continuously to the interior of the circles Si jT ei via 



u(<Si,ref) = fi. C yLt{ri,tpi) = V" ( — 



p. 

1 II 



(11) 



where r*j and are the polar coordinates with respect to circle Si >Te f and where the boundary multipoles 
satisfies P jm = -P*_ m . Note that the choice of it(5i, re f) is not unique since the continuity at the boundaries 
dSi lIe i is the only requirement. The specific choice in Eq. (|11|) is convenient for the further calculations since 
A/i )ex t = in Si^ re {\dSi^ e {. Extending the integration domain of the capillary wave Hamiltonian in Eq. 
f2 = R 2 \ UjiSj.rof — > R 2 generates an additional energy contribution 



d 2 x 



27T7 



S i>Te f 



m \Pi~ 



m> 1 



(12) 
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In Eq. I|12fl we have already omitted the contributions from the gravitational term in Ti. cw which are of order 
(r /A) 2 <C 1. From Eq. (|12fl we see that the additional terms created by the extension of the integration domain 
of the capillary wave Hamiltonian 7i cw are not constant in the - possibly fluctuating - boundary multipole 
moments Pi m and therefore lead to artificial contributions to the partition function Z. These unphysical 
contributions have to be corrected by adding TLi }COTr to the extended capillary wave Hamiltonian 7i cw [f2 = R 2 ]. 
The total Hamiltonian then reads 

2 

Wtot = W cw + 2_j [^> b + %,corr] ■ (13) 
i=l 

As in the previous sections the partition function is written as a functional integral over all possible configurations 
of the interface position u and the boundary lines, expressed by 

Z = Z^ fvuf[J Vh J[ 6[u(xi) - A«t(xi)] ™p{- nt f B T U] } ' (U) 

where the product over the S- functions enforces the pinning of the interface at the positions of the colloids. In 
contrast to Eq. l|7jl. these product extends over all x S Si,ref instead of dSi^ c {, only. The analysis of Z in this 
form for intermediate asymptotic distances d is presented below in Sec. IIII Dl 

III. LONG-RANGE BEHAVIOR 

In this section analytical expressions for the fluctuation induced force in the intermediate asymptotic regime 
r C d < A c are calculated. In Sees. 1111 Al and 1111 HI we will evaluate the partition functions of fluctuation 
and the mean- field part, respectively. From these two contributions we obtain the asymptotic form of the 
total Casimir force, which in turn will be calculated directly in Sec. IIII Dl applying the approach introduced in 

sec.nrni 

A. Fluctuation part 

The fluctuation part appears in all cases (A) and (B1)-(B3) for the boundary conditions introduced in 
Sec. Ill Al In the case (Bl) (pinned contact line with frozen colloid) it constitutes the full result for the partition 
function because u m { — and Z m { = 1. We express the 8- functions in the fluctuation part of the partition 
function l|l(J|) by their integral representation via auxiliary fields "02 ( x i) defined on the interface boundaries 
dSi iTe f. This enables us to integrate out the field u leading to 



'flue 



TT exp J V f dti [ dl, V<(xi) G(|xi - I) > , (15) 



where d£{ is the infinitesimal line segment on the circles dSi^ c {. In Eq. I|15|l we introduced the Greens function 
of the operator (— A + A^ 2 ) which is given by G(x) = ito(l x l/A c )/(27r) where Kq is the modified Bessel function 
of the second kind. In the range d/X c <§; 1 and ro/X c <C 1, we can use the asymptotic form of the Kq for 
small arguments, such that 2ttG(\x\) sa — ln(7 e |x|/ 2A C ). Here, 7 C ss 1.781972 is the Euler-Mascheroni constant 
exponentiated. 

We introduce auxiliary multipole moments as the Fourier-transforms of the auxiliary fields ipi on the contact 
line circles dSi !le {, 

$im = ^ [^ d^e^MM^)) (16) 







The analogous multipole decomposition for the Greens function G(|xj — Xj|) is calculated in App. [5] Using it, 
the double integral in the exponent of Eq. (|15fl can be written as a double sum over the Fourier components, 
consisting of a self-energy part 

<?sdf = I d£i dii ipi(xi) G(|xi - Xi|)<0i(xj) 

„■ J dSi rcf J BSi rof 
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v c 7 n>0 



(17) 



and an interaction part 

Gint = 2 



f dl x ( d* 2 Vi(xi)G(|xi -x 2 |)^ 2 (x 2 ) 

^'S'l.rcf JdS2,vcf 



2tt 



~- hi ( ^ ) ' 'Hi' '2ii 



£ 

m,n = 
n+n>l 



(-1)' 



to + ?i 
m + n \ n 



''() 



1p2n + V'l-m'02-n ) 



(18) 



Inserting these expressions into Eq. (|15fl . the functional integral over the auxiliary fields can be written as an 
integral over their multipole moments. The partition function then reads 



Z fluc = J Y[ Vipi exp • 



fc B T ( if?! 
' 2 7 



\fr. 



G S elf Gi n t 
Gi n t G sc if 




(19) 



where the vectors = (ipio, i/'iij 



.) contain the auxiliary multipole moments of colloid i. The coupling 



matrix G which contains the Fourier modes of the Greens function G(x; — Xj) has a block structure. The 
self energy submatrices G so if which describes the coupling between auxiliary moments of the same colloid are 
diagonal, and their elements can be read off Eq. I|17fl . The offdiagonal blocks Gi n t characterise the interaction 
between the multipole moments residing on different colloids and are assigned by Eq. (|18f) . Note that all matrix 
elements coupling modes with positive and negative m vanish. From Eq. (|19|) we find that the fluctuation part 
of the free energy reads 



-k^T In Z nuc oc k^T In det G 



(20) 



where in principle the matrix G is infinite dimensional and divergent and therefore requires regularisation. 
However, after performing the logarithm to calculate the free energy these divergencies only reside in terms 
independent of the geometrical arrangement of the two colloids, and therefore the derivative of the logarithm 
with respect to d which corresponds to the Casimir force is finite. Hence, the physically important properties 
of iJ nuc reside in the offdiagonal blocks of the matrix G which describe the interaction between the colloids 
depending on their separation d, whereas the self energy part ensures the correct normalisation of the interaction. 
Eq. H18|) allows for a systematic expansion of the logarithm in Eq. I|2(J|I in powers of tq / d, 



2n 



(21) 



where the coefficients f^ depend on the logarithms — ln(7 c d/ 2A C )) and — ln(7 c r / 2A C ). We show the first few 
coefficients in Tab. [I] The number of auxiliary multipoles included in the calculation of the asymptotic form of 
•7"fluc in Eq. I|21|) is determined by the desired order in ro j d. The individual contributions in this expansion can 
be understood as (possibly higher order) auxiliary multipole-multipole interactions, related to specific products 
of the matrix elements of G. So, e.g., all the logarithmic contributions to the coefficients in Eq. 1)21 [I are 
related to interactions of the auxiliary monopoles. Note that we cannot perform limA c ^oo / 2 " c because of the 
logarithmic contributions mentioned above. These logarithmic divergencies are reminiscent of the long-range 
correlations of the capillary waves which lead to a width of the free interface which diverges logarithmically 
with A c or with the system size if gravity is absent. For the Casimir force itself, however, we find a finite value 
in the limit A c — > oo. In the asymptotic range r /d <C 1 and in the limit X c /d — > oo the leading term of the 
fluctuation force is governed by the first term in the series Eq. (0, /o uc , and reads 



d 

-Ffluc = k B T— lnZfl, 
ad 

Note, however, that the limit A c — > 
l/(dlnA c ), 

d/o fluc <i/x c _«i 
dd 



fc R T 



1 



0(d 



-3\ 



d d 
- » 1, — 

ro A c 



. 



(22) 



2 dln(d/r ) 

oo here is slowly converging with a leading correction term of the order 



1 



d\n(d/ro) 



31n(r /d) 
2rn(r /A c ) 



+ O((ln(r /d) 2 /ln(r /Ac) 5 



(23) 
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TABLE I: Leading coefficients of Eq. 
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and that the free energy difference corresponding to Ea. (|22fl . T{d) ~ lnln(ro/d), is actually ill-defined, because 
the effective colloidal interaction in case (Bl) - fixed colloids and pinned interface - is only meaningful for a 
finite capillary length A c . 

B. Mean- field part 

The calculation of Z m { (Eq. (JTJJJ) requires the solution of the Euler-Lagrange equation ©, (—A + 
A~ 2 ) Umi (x) = for x G R 2 \ Uj^ref and the (fluctuating) boundary conditions 

Mmf(x 4 ) = fi(Xi) (24) 

with x, S dSi ti:c f and u m { (x) — ► for |x| — > oo. In fact a solution to a similar problem (with Neumann boundary 
conditions) was given in Ref. |38| in terms of bipolar coordinates which, however, is involved if applied to the 
general Dirichlet boundary conditions in our case. For our purpose it is more convenient to write the solution 
as a superposition 

M mf (x) = ui(x - ri) + m 2 (x - r 2 ) = ui(ri, (,ci) + u 2 (r 2 ,^ 2 ) , (25) 

where is the center position of circle Sj )re f. The general solutions itj of the mean-field equation © in polar 
coordinates (ri,cpi) with respect to the centers of the reference contact line circles Si ;Te { (see Fig. [5] ) can be 
written as linear combinations 



in 



(26) 
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The functions aj m are defined by 



Kpin/Xc) _ ln( 7e r 8 /2A c ) 
K {r /\ c ) ~ ln( 7c r /2A c ) 



OimiWi) = ##7TT ^ « M ^"^ . (27) 



#m(ro/A c ) \r, 

which are normalized for convenience and where we have used the asymptotic form of the modified Bessel 
functions K m for small arguments rj/A c <C 1. As the solution has to match the boundary conditions at both 
circles dSi lIe { and dS 2 ^ re {, we project u m { onto the complete set of functions on dSi tIe {, e tmVl , and on dS2,iel, 
e"™^ 2 , respectively. The expansion coefficients of these projections of u m { must equal the boundary multipole 
moments at the corresponding circle which leads to a system of linear equations for the expansion coefficients 
{A im ,B im } of the m, 

"-max 

Ai m + 

A 2n = Pi 

n— — n max 
n max 

A 2m + 

a 2,mn Ai n — P 2m ■ (28) 

In Eq. H28|) . the matrix coefficients ai, mn are given by the projection coefficients 



1 f 2 * 

oi, m » = ^ J d<pie mifil a 2n {r 2 ,^ 2 ) , (29) 

and analogously a 2 , mn is given in terms of the projection of ai m (ri,ipi) onto e mlip2 . In Eq. I|28|l. we have 
truncated the Fourier expansions at a maximum order n max - Indeed, the numerical solution shows rapid 
convergence of the expansions in Eq. (|26[) in the whole range of colloid-colloid separations d, such that they can 
be truncated at n max ~ 20. Through the inversion of the linear system (|28|l the coefficients Aj m and Bi m can 
be expressed as linear combinations of the boundary multipoles, 

2 n m „ 

Aim — ^ y ^ ^ Pimjn-Pjn ■ (^0) 
j=l Tt=—rttxxax 

Using the asymptotic forms of the functions ai m given in Eq. I|27|l . expressing r 2 — r 2 (d,ri,ipi) and (p 2 = 
l p 2 {d, ri, (pi) by the polar coordinates with respect to circle Si. rc f and vice versa for circle <S2,refj an d expanding 
the projection coefficients ai :inn in Eq. (|28|l in ro/d, the coefficients Pi m jn in Eq. 1|30[1 can be written as a power 
series in ro/d. Inserting the solution (|25fl with the expression 1|30|) for the Ai m into the line integrals for the 
mean field energy 7Y[it m f] in the exponent of Eq. I|10|l . we can write 7i[u m f] as a quadratic form of the boundary 
multipole moments Pi m , described by a matrix E with block structure, 




(31) 

Here, the submatrices E,, se if and E; n t describe the coupling energy of the multipole moments of the same contact 
line fi and from different contact lines, respectivly, and their elements are given by 



E TZ\i = 27r5 rn „(l - S m o) - / rod^e"^ 1 2,Pikin-^-a ik {ri,ipi) , (32) 

JO k 0r i 



£/ 12.int 



/•2-rr Q 

/ r dipi e lmVl y^P2kin-^—a2k{r 2 ((fii),ip 2 {(pi) , (33) 



and analogously for S^i'int ■ Applying again the expansions of the functions a,-& in powers of ro/d in the expression 
for E^SX-at an( i using the series expressions for the coefficients Pi m jn, the matrix elements can be written as a 
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TABLE II: Leading coefficients 2/ 2 ™ f /k B T of Eq. JH9 for various boundary conditions. The coefficients /^ uc are given 
in Tab-E 



2n 


case (B2): only monopoles 


case (B3): mono- and dipoles 


case (A): all multipoles 
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/"flue 
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/"flue 
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/"flue 
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-f2 uc - 2 


/"flue 
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/"flue 
— ./4 
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-/ 6 fluc - 8 


/■flue 
— JO 


/"flue 
— ./6 


8 


-fi uc - 31 


-fi uc - 18 


-/I" - 2 



power series in ro/d in the asymptotic range r$ <C d <C A c . The mean field part of the partition sum Eq. (|10|) 
can then be written as 

such that the d-dependent part of Z m { is proportional to detE. As discussed in Sec. Ill Al the functional 
measure is given by a product T>fi = dhiY\^__ M dPi, n , where M depends on the choosen model for the 
boundary condition at the contact line. Inserting the power series of the matrix elements E™™ t and -E™"i f and 
expanding the free energy J- m f(d) oc /cjsTlndetE, we arrive at a series similar in form to the fluctuation part, 
Eq. :EB- 



n 

where again the coefficients are functions of the monopole self-energy and interaction coupling elements, 

— ln(7 c r /2A c ) and — ln(7 e d/2A c ), respectively. The exact form of the can be recovered from Tab. [D] Before 
discussing the analytic structure of the dependence on d, we remind the reader that the integration measure 
Vfi differs between the cases (A) - no pinning, (B2) - pinning and height fluctuation of the colloids, and (B3) 

- pinning with collective height and tilt fluctuations of the contact line. 

In all these cases for the boundary conditions we have a mean-field contribution, and the leading term of the 
free energy in the long-range regime d 3> ro (with A c — > oo) is determined by the boundary monopole-monopole 
interaction and leads to a repulsive effective force between the colloids: 

„ , „ d knT 1 o, d d , . 

F mi = k B T— lnZ mf -> | + Of 3 , -» 1.-T--0. 36 

od 2 d\n{d/r ) r A c 

So the leading term of the mean-field contribution is equal in size but opposite in sign to the fluctuation 
contributions, and thus these terms cancel. We will discuss the leading terms of the total Casimir and their 
connection to the type of boundary conditions which is considered in the next subsection in more detail. 



C. Total Casimir force for d ^> ro 



The 1/d asymptotics for both the mean-field and the fluctuation part at hand, we will discuss the total 
effective force F, which due to Z = Zf[ uc Z m { is obtained as F = F fluc + F mf . As discussed before, the mean- field 
result depends on the type of boundary conditions under consideration. In all cases the large d expansion of 
the mean-field part of the partition sum (with A c — ► oo) leads to a repulsive effective force between the colloids, 
whereas the fluctuation contribution is attractive. So the total Casimir force between the colloids is determined 
by an interesting interplay between the two-dimensional "bulk" fluctuations and one-dimensional boundary 
fluctuations influencing the two-dimensional "bulk" by a change of the mean-field. 

Combining the expansions of the mean-field and the fluctuation free energies, Eqs. I|35|l and l|21|) . the asymp- 
totic form of the total Casimir can be written as a power series in ro / d, 

™--¥!£[s?+/Sfl(?f ■ m 
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For a pinned contact line and fixed colloids, = for all n, and the full Casimir force is given by the 
fluctuation part result Eq. (|22p. As we have seen above, for fluctuating boundary conditions, the leading terms 
of F m f and Fa uc (Eqs. (|36(1 and l|22|l) cancel. The same holds for the first subleading terms, for which we find 
uc = — /™ f for all types of boundary conditions with a mean-field contribution. The form of /™ f , however, 
is different for the cases (B2), and (B3) and (A), respectively. In the simplest case of a pinned contact line 
with height (or monopole) fluctuations of the colloid, we find /| nf + ff uc = —1, such that in case (Bl) the 

leading term of the fluctuation-induced force in the limit ^ > 1 , j > is attractive and stemming from the 

fluctuation part, and can written as 

F^-— 5- (4) . (38) 



r \d 

For a pinned contact line with tilt and height fluctuations (boundary mono- and dipoles, case (B3)) and for 
the generic case of an unpinned and fluctuating contact line (described by all boundary multipoles, case (A)), 
we find that this term also cancels: f™ 1 = —f^ and /™ f = — /g . In these cases, the first non-vanishing 
coefficient in the expansion Eq. (|37|) is stemming from the fluctuation part as well. In case (B3) it is given by 
/s mf + fi uc = - 9 which leads to 



For case (A) we find /™ f + /f uc = — 1 leading to an attractive Casimir force which (in the limit ^ > 1 , j ► 0) 

reads 



D. Kardar's method 



Here we calculate the partition function directly without splitting the fluctuating field u into a mean-field 
and a fluctuation part. The starting point is Eq. Ijl4|) from Sec. Ill CI The (5-functions can again be expressed 
by auxiliary fields tpi, now defined on the two-dimensional circular domains S^ )re f as opposed to the auxiliary 
fields of Sec. IIIIA1 which are defined on the one-dimensional circles dSi xe f- 



2 

Z = [Vu f[V^ Vfi exp 



Htot[/i,M] 



ii d xipi(y:)[u(x) - /j, ex t(x)] 

'Si.ref 



(41) 



The function /^ext which describe the interface extension to iSj jre f have been defined in Eq. (|ll|l in terms of a 
multipolc expansion in the Pj m . The total Hamiltonian contains the capillary wave Hamiltonian, the boundary 
and correction terms (see Sec. Ill C|) and reads 



Htot = I I d 2 x 



(Vu) 



A? 



TT", 

T 



E 



|m|>l 



(42) 



Similarly to the evaluation of the fluctuation part, Sec. IIII Al we introduce multipole moments of the 
auxiliary fields by inserting unity into iJ, Eq. (|41|l : 



2 

J i=l m 



d 2 a;(r/ro) |m| e- imv V(x) 



(43) 



In contrast to the evaluation of the fluctuation term in Sec. 1111 Al there will be constraints on the lowest 
multipoles which contribute to Z. To see this we note that, after shifting hi — > hi — Pm, the Hamiltonian Titot 
does not depend anymore on the boundary monopole moments Pto and the dipole moments Pn and the only 
dependence of Z on these moments is through the constraint function fi iCX t- Recalling the definition of the 
integration measure D/j for the various boundary conditions, Sec. Ill Al and performing the integration over P i0 
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and Pn where applicable, we immediately find 

2 



(44) 



Jf[Y[d9 im ...5(9io)... case(B2) 

2=1 m 

r 2 

/ II II • • ■ S ^^) S(*i-i) • ■ • cases (A) and (B3) 

•* i=l m 

Having noticed these constraints on the auxiliary fields, we proceed by integrating over the field u in Eq. I|41|) : 
Z = I JJl>^i / Vfj cxp < ^ / d 2 ^i / d 2 ^- ^(xi) G(|xi - Xj-|) ^-(Xj) 

~*E^ d 2 x^(x)/., oxt (x) 1 , (45) 



7T7 



2(P i0 -^) 2 + 4 J] (1- H)|P„ 



where - as in Eq (|15fl - G is the Greens function of the capillary wave Hamiltonian. A somewhat longer 
calculation shows that Z can be split into into an interaction part (coupling the auxiliary multipolc moments 
^>i m for different colloid labels i), a self-energy part (depending on for each value of i separately) and a 
remainder (the sum of boundary and correction Hamiltonian): 

Z = j nn^™ / V ^ eXP {"T^ ( W M[ $ l""*2m] +Hi, ae u[% m ]) 
i=l m •* I - 



exp 



7T7 



2(P M -h l f+ (l-H)|Pi, 

|m|>l 



(46) 



The interaction part 



2 j d 2 Xl f <Px 2 Mxi)G(\xi -x 2 |)^ 2 (x 2 ) 

JSl >te { J S2,rof 



1 

2^ 



-2 In 



7c d 
2\ r 



10*20 



E 

l± ,l2=° 
l!+l 2 >0 



(-1)' 2 fh + h\ fr 



h+h\ I 



/j"o\ il +' 2 

\~d) (*l'i* 2 'a *2-( 2 ) 



(47) 



is a bilinear form in the auxiliary multipolc moments; it was derived using the multipole expansion of the Greens 
function G(|xi — x 2 |) ~ — ln(7 e |xi — x 2 |/2A c ) (valid for d 3> ro) which is presented in Add. 151 in more detail. 
The self-energy part 



exp 



k B T 
' 2 7 



TU^oj = J II ^ ^ (*im ~ f s d 2 x (r/r ) |m| e- im ^V(x)) exp ^ £ im P it 



(48) 



exp 



k B T 

' 2 7 y Si 



cf a; / d 2 x't(j l (x)G(\x-x.'\)t(j i (x) -i dfx ^i(x) /i, ext (x) 



is evaluated in App. O with the result 



Si 



2 ln( 7e ro/2A c ) ^ |J 



2tt 



E 

m>0 



27r|m| 



Combining Eqs. I|46|) . i|47|l . and Ij49(l . the partition function can be written as 



Haelf Hi nt 

Hgclf 



*1 

* 2 



(49) 



(50) 
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where the vectors *S?i — (^io, Pm, ^ii,Pn, Pi-i, . . .) - in contrast to \&i in Sec. lIII Al - contain all involved 
auxiliary and boundary multipole moments. The elements of the matrix H describe the coupling of these 
multipole moments, where the self-energy block couples multipoles defined on the same circles S^ref- The 
self energy matrix H so if can be read off Eqs. ((46(1 and 1)49(1. The elements of the interaction matrix Hj llt are 
determined by the interaction energy 7i m t in Eq. (|47|l and couple the auxiliary multipole moments of different 
colloids. All matrix elements representing couplings of other multipoles are zero. 

Similar as in Eqs. 1(19(1 and 134(1 . the exponent in Eq. ((50(1 is a bilinear form, however, here combined for 
all types, boundary multipole moments Pi m and auxiliary multipoles ^i m . The computation of the partition 
function amounts to the calculation of detH. Expanding the logarithm of this determinant for r^/d <C 1, and 
taking the derivative with respect to d, we directly obtain the asymptotic form Eq. 1(37(1 for the total Casimir 
force, 

n 

with /2« rdar — fin C + fin as ^ should be. In contrast to the calculation before, the different leading power laws 
for the different cases (A) and (B1)-(B3) can be understood easily. We note that the interaction between the 
auxiliary multipoles fi m and in Ti- m t, Eq. ((47|) . scales like (tq / d) m+n . Thus, the leading order of the total 
fluctuation induced force between the two colloids is determined by the first non-vanishing auxiliary multipole 
moment ^i m ' and (as follows from det H) gives rise to a force F(d) oc l/d 2m +1 (for to' > 0) or F(d) oc l/(dlnd) 
(for to' = 0) . As explained in the beginning of this subsection, the different boundary conditions lead to certain 
constraints on the auxiliary multipoles: According to Eq. 1(44(1 . the leading term in F(d) arises from a monopole- 
monopole interaction of the auxiliary field in case (Bl), from a dipole-dipole interaction in case (B2), and from a 
quadrupole-quadrupole interaction in cases (B3) and (A). We remind the reader that the constraints of vanishing 
auxiliary monopole and dipole moments result from the independence of Titot of the boundary monopole and 
dipole moments and that this is only captured correctly by the inclusion of the correction Hamiltonian 7i corr 
(see Sec. HTTP)! . 

At this point we insert the following observation: If the sum of boundary and correction Hamiltonian were 
zero, all multipole moments f ; ra would be zero and consequently all coefficients in the expansion of the Casimir 
force in terms of r^/d would vanish - i.e. the total Casimir force would be zero. This happens for the boundary 
Hamiltonian [3jJ 

W b ,i = 2^ 7 ^m|P„„| 2 (52) 

m>l 

7 r ^ r ,m m^-hm? (53) 



l ^J dSi vct J aSt ^ sin 2 [\ (<£-<£')] 

Thus we see that the boundary Hamiltonian needs to be of nonlocal nature in the contact line height fi to make 
the Casimir force vanish. 



IV. SHORT-RANGE BEHAVIOR 
A. Fluctuation part 

In the opposite limit of small surface-to-surface distance h = d — 2r$ <C tq the fluctuation force can be 
calculated by using the Derjaguin (or proximity) approximation |3lj . It consists in replacing the local force 
density on the contact lines by the result for the fluctuation force per length f2d(h) between two parallel lines 
with a separation distance h and integrating over the two opposite contact line half-circles to obtain the total 
effective force between the colloids. 

The Casimir force between two parallel surfaces was calculated in Ref. |3r| in a general approach and explicitly 
for three-dimensional problems. Applied to two dimensions we obtain f2d(h) — — k-gT tt/ (2Ah 2 ). Integrating 
over the opposing contact line half-circles yields 

77 ~ O^B? 1 f r ° , 1 ro/h^oo 7T 2 rl /2 nfi-Xlls /~,N 

-Ffluc « -2^j— / dy ; g — > -k B T — + 0{h 1 ) . (54) 

24 Jo (h + 2r -2^W) 48/13/2 
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for the dominant contribution to the Casimir force from the fluctuation part in the limit h/ro <C 1- Note, that 
this strong increase as h — > is a consequence of the finite (mesoscopic) size of the colloids and is missed if the 
colloids are approximated by pointlike objects |4flj . 



B. Mean- field part 

Here, we exemplify the asymptotic behaviour of the mean-field force for close colloid separations h by 
case (B2) for the boundary conditions where only fluctuations of the boundary monopoles occur. From the 
numerical results (see next section) we observe that the type of divergence of the mean-field force as h — > is 
obtained correctly by considering only monopole fluctuations; including higher multipole moments affects the 
force only by a multiplicative constant. (This is in marked contrast to the long-range regime.) 

In order to apply the Derjaguin approximation, we calculate the mean- field between two parallel lines dSi/ 2 
on which the field is pinned to a fluctuating value u m {(dSi) — P,o (corresponding to monopole boundary 
conditions). The mean-field energy 7i m i in Eq. <|31[1 is represented by a 2 x 2-matrix E. By diagonalizing 
E we can write H m f — ^-^-{e^Z 2 + e a (Az 2 )), where e s and e a are the eigenvalues of E, and Z = P w + P 2 o 
and Az = P w — P 2 o are the symmetric and the antisymmetric superposition of u mf at the boundaries dSi, 
respectively. L y is the length of the boundary lines. In fact, the line densities e s and e a correspond up to a 
factor to the mean field energies of the solutions of the mean-field equation (|24|l with symmetric (w s |<9 s = Z) 

1/2 

and antisymmetric (u a |9 s = ±Az) boundary conditions, respectively, and the general mean field solution is 

1/2 

given by by u m f — (u s + u a )/2. For A c — > oo the mean-field equation reduces to the Laplace equation Au m { = 
between parallel lines which are in x-direction a distance h apart. The antisymmetric and symmetric solutions 
read w a = (2Az/h)x and u s = const., respectively, with the corresponding line densities e s — and e a = i/h. 
A vanishing e s signifies that a collective vertical shift of the interface does not cost any energy. This leads to the 
divergence of the integral over Z in the partition function Z m {, which physically is not interesting and can be 
neglected. Applying the Derjaguin approximation (similar to Eq. 124(1 ') to only the antisymmetric mode yields 
the corresponding energy E a for the two circles a distance h apart: 

E a « -2 7 P dy , ^ 2(A f ] ] Q o re 27r 7 (Az)\f^ + 0(1) . (55) 



h + 2r - 2vVg - y 2 ' V h 



Thus the h dependent mean-field part of the free energy reads J-'uif — —kBTln J d(Az) exp(— E a j k^T) = 
(fceT/4) \n(h/r ) + const. So, in the limit h — d — 2R — > the leading (divergent) part of the effective mean- 
field force F m { is repulsive and reads 

F m{ (h^0)^^. (56) 

It appears physically less obvious why the Derjaguin approximation could also be applied to higher boundary 
multipole moments n, especially if h < R/n. If one nevertheless does so one finds that in this regime h < R/n 
the mean-field force diverges slower than 1/h whereas for h -C R/n the monopole behavior is recovered. This 
is in accordance with our numerical results. 



V. INTERMEDIATE DISTANCES: NUMERICAL CALCULATION 



For intermediate distances d ~ ro the fluctuation induced force has to be calculated numerically We will 
do this as in the previous sections for the fluctuation and the mean-field part separately. For the fluctuation 
part, we shall apply a method which was introduced in Ref. [l^- The starting point is Eq. (|15|) for the 
partition function Zfl uc . Introducing an equidistant mesh with TV points <pij = (2n/N)j, < j < N, on the 
contact line circles dSi iie { converts the double integral in the exponent to a double sum. Then the functional 
integrals over the ipi are replaced by ordinary Gaussian integrals over the tpi(xi((fij)) , T>ipi ~ YijLo dipi{xi{<fij)) ■ 

In the exponent, the ^(x^^-)) are coupled by a matrix G with elements G J J, = G(|xj(<^j_y) — x$/(^/,-/)|). 
Performing the Gaussian integrals and disregarding divergent and (i-independent terms immediately leads to 
fj uc = (k?,T/2) lndet(G^ c 1 G(d)) for the fluctuation free energy. Here, Goo = limd->oo G(d). It contains the 
self energy contributions and is needed for the regularization of the free energy. Deriving with respect to d, the 
Casimir force can be written as 

F nuc (d) = -^tv[G(d)- 1 d d G(d)] . (57) 



1G 



The advantage of the direct calculation of the force is that Eq. I|57|l does not contain any divergent parts which 
would require regularization, thus easin g th e numerical treatment considerably. The determinant is computed 
by using a standard LU decomposition jij. We find good convergence of the numerical routine also for small 
d if N « 5000, which, however, demands a computation time of about 30 hours for each distance point on a 
standard PC. As discussed in Sec. I VII and shown in Figs. we find very good agreement between analytical 
and numerical results. 

The numerical calculation of the mean-field Casimir force can be done very conveniently with the method 
described in Sec. lIII 61 In order to avoid complex numbers, we used sin- and cos-modes instead of Fourier-modes 
for the numerical computations. It is straightforward to rewrite the corresponding equations in Sec. IIII Bl by 
using their real and imaginary parts. Via a numerical evaluation of the integrals for the projection coefficients 
ai i7nn in Eq. I)29|) and inversion of the linear system Eq. (|28|) the elements E™£ 1{ and E^ n int of the matrix 
E are computed, see Eqs. Ip?2")) and (|3*5|) . Calculating IndetE then provides the free energy ,F m f, and the 
numerical derivative finally the mean-field Casimir force F m [. This scheme turns out to be very efficient and 
provides results for F m [ within seconds. Indeed, the numerical calculations show that it is sufficient to consider 
n mm ss 20 modes in the expansions of the Ui in Eq. i|25|) in order to achieve convergence of the results. 

VI. DISCUSSION OF RESULTS AND OF POSSIBLE REALIZATIONS 

In Fig. we compare numerical results with the analytical expressions for the Casimir force for the long-range 
asymptotics fSec. IIIlH)) and for very small colloid-colloid separations fSec. IIV(I . respectively. 

For all cases considered for the boundary conditions, (B1)-(B3) and (A), we find a very good agreement with 
the analytical predictions, Eqs. (|22|l and l|38() - (|40|l in the long range regime for colloid separations d > 5ro. 
In the short-range regime, the fluctuation induced force depends much less on the specific type of boundary 
conditions at the contact line. In fact, in this regime, the divergence of fluctuation part force (which applies 
to all cases), ~ —h~ 3 / 2 , is dominating the total force in this regime and leads to a strong net attraction. This 
may have an important influence on the effect of colloid aggregation at interfaces as is discussed in more detail 
in Ref. [2l| , especially with regard to the van der Waals force. As can be seen from the plots in Fig. [3] in case 
(Bl) (fixed colloids, no mean-field force) the Derjaguin approximation describes well the numerical data also 
for intermediate colloid separations r$/h = ro/(d — 2ro) < 5. In the plots for cases (B2), (B3) and (A), however, 
we see that the repulsive contribution from the mean-field part showing a weaker divergence ~ 1/h leads to a 
strong decrease of the total force for r^/h = ?*o/(d — 2ro) > 1, which, by increasing d, converges to the power 
laws governing the long-range asymptotics. The rapid decrease of the total force actually renders the Casimir 
force effectively short-ranged for the cases (B3) and (A). Note that because of this rapid decay ~ —1/d 9 , the 
competing attractive and repulsive contributions from the fluctuation and mean-field part are of almost equal 
size and our numerical methods which are based on the addition of these quantities are afflicted with numerical 
uncertanties and are not able to provide reliable results for d > 15ro. 

In experimental realizations, fixing the colloids as required in case (Bl) might, e.g., be realized by a laser 
tweezer. In such a setup, the fixing is usually, of course, not exact, but means that the vertical movement of 
the colloid centers is restricted by an external potential. Including such a external potential for the colloids in 
our model in fact lowers the repulsive contribution from the mean- field part (as compared to the unfixed case) , 
and therefore leads to an increased long-ranged attractive total Casimir force dominated by the leading term of 
the fluctuation part, Eq. (|22|l . The extent of this increase is controlled by the strength of the external potential 
for the colloids, as we will show in the next subsection. 

A. External potential on the colloids 

In the effective Hamiltonian introduced in Sec. [n] we considered only free energy differences resulting from 
the changes of interfacial areas which are associated with fluctuations around the flat reference configuration. 
In this subsection we will extend this model to external potentials Vi(hi) acting on the vertical position hi of the 
center of colloid i. We will concentrate on the cases of a constant external force Fie z and a harmonic potential 
for colloid i, corresponding to Vi(hi) = —Fihi and Vi(hi) = ^Di (hi — h ^) 2 , respectively. The vertical forces Fi 
include, e.g., gravity, whereas the harmonic potential can be realized by a laser tweezer. So this is of particular 
interest for case (Bl) where the colloid positions are fixed. 





d/r d/r 

FIG. 3: Comparison of the numerical results for the Casimir force (symbols) with the analytical expressions in the 
asymptotic ranges of large colloid separations d 3> ro (full line) and small surface-to-surface distance h — d — 2ro <C ro 
(dashed line) for the boundary condition cases (Bl), (B2), (B3), and (A), respectively. 



Then, the total effective Hamiltonian of the two colloids adsorbed at the fluid interface reads 

n = n CVJ + J2 {«b,i + v t (ht)} . (58) 

i 

Note, that the approximations performed in Sec. [n] in deriving the capillary wave and boundary Hamiltonian 
remain valid here if the external force or the displacement of the harmonic potential are small on the scale set 
by the surface tension, i.e., Fj <C and Dih <C 27T7r , respectively. 

The additional external potential has two implications which we will discuss in the following. First, it leads 
to a deformed equilibrium meniscus as compared to the flat reference interface, which gives rise to a "classical" 
capillary interaction between the colloids [I3 ■ Secondly, through the coupling of the colloid position hi to the 
interface field u in the boundary Hamiltonian 7ib,i, the thermal movement of the colloids in the potential Vi 
can also influence the fluctuation induced ("non classical") force between the colloids. 

As described in detail in Ref |l3j . the equilibrium meniscus u oq can be found by minimizing the effective 
Hamiltonian H with respect to the colloid position hi and the interface height u(x,y). The equilibrium colloid 
height ft-^cq is determined from the condition dTl/dhi = and depends on both the external potential and 
the mean height of the three phase contact line on the colloid surfaces. The interface field u cq fulfills the 
Euler-Lagrange equation (|24|l of H cw with the boundary condition [l3| 

du cq (yi) tt eq (x) - fe Mq 
dru r 



(59) 



18 



where d/drii is the derivative in the outward normal direction of dSi^rf. Using the general form of Eqs. (|25|) 
and l|2rj)l for u eq and projecting the boundary conditions on dS^rcf onto e lmipi as in Sec. 1111 Ml leads again 
to a linear system for the expansion coefficients Ai m similar to Eq. (|28|1 . 

For a constant external force, we obtain in the asymptotic range rg <C d <C A c 

Ucq ^ _Ai n f^ _ I^o ( r_o\ cog( } _ ^L ln _ (r_o\ cog( } (flQ) 

q 2tt 7 V 2A c/ 2vT7dV? , i/ 2tt7 \ 2A C / 2wjd \r 2 J 

for the leading terms in ro/d of the equilibrium meniscus u oq (see Fig. [21 for the definitions of r t and ipi). The 
capillary interaction arising from the meniscus deformation is given by 



V men (d) = H(d)-J2 H i 



2irj 2A C 47T7 



Here, Tii :OQ is the effective energy associated with a single colloid system. The first term in the second line 
of Eq. IjtjlD is the well known logarithmic flotation force. Note, that in the absence of electrostatic forces the 
main contribution to the external force usually is gravity which can be neglected for colloid radii < 1/im as it is 
much smaller than the thermal energy k^T. The second term describes the leading behaviour of the capillary 
interaction if there is an external force only applied to one of the colloids - and it leads to a repulsive force 
between the colloids. 

For harmonic external potentials the equilibrium meniscus and the capillary interaction potential is calculated 
in the same way but results in a somewhat lengthy expression for V mon . Focusing onto the symmetric case of 
identical potentials for both colloids, we find that Vmon vanishes in the limit A c — > oo: In the absence of gravity 
a parallel shift of the whole interface does not cost any energy, and, therefore, the response of the system on the 
presence of the harmonic potentials on the colloids is a shift of the planar interface to a new equilibrium position 
u eq = — ho (= hoi) - without a meniscus deformation and, hence, without causing a capillary interaction. The 
leading term for 1/A C small is given by 

v 4ir 7 fegln( 7c d/2A c ) 

Vmcn ~ (2 7 r7/ J D-ln(7ero/2A c ))(27r 7 / J D-ln( 7e rod/4A2)) ' { ° ] 

If the minimum of the harmonic potentials on the colloids are not identical, we find an additional term 



a c ^oo 7r 7 (/ioi - h< 



■tr /v c — ,, , v .m '"02/ Ip.'W 

mcn * l + Tr^l/A + l/AO + lnd/ro ( j 

which is not vanishing in the limit A c — > oo. 

The influence of the external potential on the fluctuation induced force is most conveniently discussed by 
splitting the interfacial field into an equilibrium and a fluctuation part, u = u cq + Ait, similar to the splitting 
into a mean-field part u m f and a fluctuation part v described in Sec. Ill Bl and analogously hi — /ij ieq + 
Ahi for the vertical position of the colloid centers. Note, however, that u eq fulfills the equilibrium boundary 
conditions Q59J1. whereas u m i is the mean- field part of Au, i.e. Au = u m { + v, which fullfills the boundary 
conditions Eq. @ corresponding to thermal fluctuations of the contact line around its equilibrium position. 
Inserting this decomposition of the the interfacial field u into the extended Hamiltonian Eq. (|58|l and performing 
some conversions of the integrals by exploiting Eqs. I)24|) and (|59|l together with Gauss' theorem, we find a quite 
distinct behaviour for the cases of a constant external force and of a harmonic potential. In the first case we 
can rewrite the effective Hamiltonian in the form H[u cq + Au] = H[u cq ] +H cw [Au] +Hb[Au, Ahi]. That means 
that the effective Hamiltonian relevant for the fluctuation induced force is independent of the external force in 
this case. The only effect of a constant external force, hence, is the deformation of the equilibrium meniscus 
which leads to the classical interaction described by Eq. I|61(l . In the case of a harmonic external potential, 
however, we find H[u cq + Au] = H[u cq ] + H C w[Au] + Hb[Au, Ahi] + ^Z i {F>i/2){Ahi) 2 , i.e. compared to the 
functional integral Eq. (JJJ we have to introduce an additional harmonic potential term for the deviation Ahi of 
the colloid position from its equilibrium value /i^ e q in the expression for the partition function describing the 
thermal fluctuations of the interface and the colloids. The additional potential term Vi = (Di/2)(Ahi) 2 , which 
is centered at Ahi — 0, can be included in the boundary Hamiltonians 7i^b and leads to additional terms in the 
self-energy parts of Eq. 1-iL'l) for the mean-field part of the partition function. The integrals over the Ahi and 
the boundary multipoles can be performed as before. The additional potentials Vi leads to modifications of the 
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determinant of the matrix E. Because of the different form of integration measure for the cases of a pinned und 
an unpinned contact line (c.f. Sec. Ill A|l . these modifations differ for the two types of boundary conditions. For 
the leading term of the mean-field Casimir force we can write 

a cZ? oo 1 (f , A s 

mf 2 rflnrf/r 0l , o lO2 ln(rf/r ) ' [GV 

where Oi = Di/(ir"{) for a pinned contact line and Oi = 2/(1 + 2-K^jDi) for an unpinned contact line. So the 
inclusion of the harmonic potential leads to a decreased mean-field contribution to the fluctuation induced force 
as compared to Eq. (|36)1 which decreases with increasing strength Di of the potentials. For a pinned contact line, 
in this way the leading term of mean-field contribution to the Casimir force can be switched on in a controlled 
way by an external laser tweezer potential. In the limit Di — > we recover for both cases the result from 
Sec. IIII 51 Eq. (|36|l . which gives rise to the cancellation of the leading terms with that of the fluctuation part. 
This is expected because Di = corresponds to the cases (A), (B2) or (B3) where height fluctuations of the 
colloids are not suppressed by external potentials. In the opposite limit Di — > oo, however, we find F m { — > 
from Eq. H64|) for a pinned contact line. This corresponds to case (Bl) of frozen colloid positions and a pinned 
contact line and where the full Casimir force is given by the fluctuation part Eq. 1)221) . In the case of an unpinned 
contact line, this effect is less pronounced. Then the original form of the mean- field Casimir force, Eq. (|3fc>|) . is 
diminished by an additional factor 1/(1 + Ind/ro) in the limit Di — » oo. 

For the linear external potential (constant force on the colloids) - also relevant for a laser tweezer with 
relatively large displacement - the only effect is an additional capillary interaction to the fluctuation induced 
force. The relative importance of these two contributions to the colloid interaction is determined by their 
amplitudes Ff/'y anc ^ ^bT, respectively. For a harmonic external potential, the relevance of the external 
potential on both the classical capillary interaction (cf. Eqs. (|61[) -ilti3 [) ) and on the fluctuation induced force 
Eq. I|64l) , is governed by the the ratio 'y/Di of the surface tension and the stiffness , Di, of the harmonic potential. 
For typical values for the surface tension in the range of 7 ~ 10 mN/m and a stiffness Di ~ 10~ 3 — 1 mN/m 
|42j of the harmonic tweezer potential, we find that the classical interaction can be neglected compared to the 
fluctuation force (whose strength is governed by k&T), but, also, the effect of the external potential on the 
Casimir force is marginal. For steeper tweezer potentials or lower surface tensions (c.f., e.g., Ref. 0]), sucn 
that "/ / Di < 1, we find an increasing influence of the the external potential. Then the mean-field part f m f 
of the fluctuation induced force is considerably decreased which leads to a more attractive total Casimir force 
because the fluctuation part then will dominate F m f (c.f. Sec. IIII C|) . On the other hand, in this regime also the 
classical capillary interaction reaches the magnitude of the fluctuation induced force. 



VII. SUMMARY 



The restrictions that two rotational symmetric colloids trapped at a fluid interface impose on the thermally 
excited interfacial fluctuations (capillary waves) by their sheer presence lead to a thermal Casimir force. This 
effective fluctuation force can be calculated by a decomposition of the interfacial partition function into a mean- 
field and a fluctuation part both numerically in the whole range of colloid separations d and analytically for 
the asymptotic ranges of either small or large separations. In the long-range limit, d 3> ro, the resulting force 
depends crucially on the boundary conditions at the three-phase contact line on the colloids' surfaces. There one 
observes an interesting interplay between the attractive interaction from the interface fluctuations and a repulsive 
interaction caused by the fluctuating boundary conditions (mean- field part). This results in a cancellation of the 
leading terms up to a certain order in \/d, which is determined by the specific model considered for the boundary 
conditions. For freely fluctuating colloids either with a pinned (B3) or unpinned (A) contact line - the Casimir 
force is characterized by a fast decay oc — d~ 9 . For a pinned contact line, fixing colloidal degrees of freedom 
leads to longer-ranged forces which are oc — d~ 5 if the orientation of the colloids is fixed and oc — l/d\nd if both 
orientation and vertical position are fixed. This cancellation of the leading terms from the mean-field and the 
fluctuation part can be understood in an alternative approach to compute the partition function of the system 
in which the analogy of the effective Hamiltonian to electrostatics is exploited. In this approach the fluctuation 
induced force can be interpreted in terms of an interaction between auxiliary multipolc moments defined on the 
area enclosed by the contact line of the colloids. The various boundary conditions at the three phase contact 
line translate into conditions for the auxiliary multipoles. The asymptotics of the fluctuation induced force for 
intermediate colloid separations then is determined by the interaction of the leading non-vanishing auxiliary 
multipoles. 
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FIG. 4: Numerical calculation of mean-field part F m { of the Casimir force with harmonic external potential for the 
colloids for Di/j — 1 (triangles, compared to the analytical result Eq. (16411 solid line) and -Di/7 = 0.1 (pluses) in the 
case of pinned contact line. The circles show the fluctuation part — Fa uc of the Casimir force, which has to be added to 
F m f to obtain the total force. As can be seen from the plot, for increasing -Di/7 the fluctuation part becomes dominant. 

In the opposite limit of a close colloid-colloid separation h = d — 2tq <C ro, the effect of the boundary 
conditions is much less pronounced. Both the mean-field and the fluctuation part diverge for h — + 0, but the 
resulting force is dominated by Ffl uc ~ /i~ 3 / 2 (compared to F^f ~ h~ 1 ), leading to a strong Casimir interaction 
in this regime (c.f. Fig. [3J. 

For typical values in experimental situations, the fluctuation force will dominate "classical" capillary forces 
arising from meniscus deformations by an external potential. Nevertheless, such external potentials provide 
the possibility to tune the fluctuation force directly or to superimpose "classical" and fluctuation forces by a 
sophisticated choice of the experimental setting. 

APPENDIX A: THE BOUNDARY HAMILTONIAN H i b 

From Eq. |0J we find that the boundary Hamiltonians Tii.h consists of three parts. In this appendix we 
show how Hij, can be expressed in terms of the boundary multipole moments Pi m , leading to the second order 
expansion given in Eq. (0J. For the case (A) of spherical colloids and a fluctuating contact line all areas AAi/u 
and AA pro j are nonzero. For the case (B3) of a pinned contact line (disks or Janus spheres), AAi/u = and 
the boundary Hamiltonians are fully determined by the change in the projected meniscus area A74 pro j. (For 
the remaining cases (Bl) and (B2) all area changes vanish and thus the boundary Hamiltonian is zero.) It is 
sufficient to determine the area changes AA' I ^ U — and A^4p ro j for a single colloid; the total area changes are 
just given by a sum of these. 

If the three phase contact is slowly varying without overhangs it can be written as a function of polar angle 
tp, and for sperical colloids its distance to the z-axis is given by 

ro(<p) = y/R? - (ho +u(r (<p)) 2 = ^/rg - 2 R cos 9 u(r (<p)) - u(r {<p)) 2 , (Al) 

where u(ro(y>)) is the height of the three phase contact line, and ho — —Rcos8 is the height of the colloid center 
in the reference configurationjc.f. Sec.^). For the second equality we have have used ro = Rs'm9. 

Following App. A in Ref. [T3, we parametrize the projection of the actual three phase contact line onto the 
reference plane in terms of the polar angle tp and write {h is the change in vertical position of the colloid center 
with respect to the reference configuration) 

7iA^ + 'YnAA'n - 1 jj d v [u(r (<p)) -h] 2 + l jj d V [r 2 (^) - r 2 ] 
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- 2 fj ^ U ~ k]2+ 2 C ^ " ' (A2) 

where in the first term in second line we have replaced u(ro(ip)) ~ u{tq, (ft) = /, i.e. we have replaced the contact 
line height by the meniscus height at the reference contact circle, since this approximation produces only terms 
which are at most of third order in the boundary multipoles. The second contribution to Tib stems from the 
changes in the projected meniscus area and can be written as 

,2-7T /TO /-27T 

7 AA' proj - 7 / dtp drr = l dp [rg - rjjfc)] (A3) 

JO Jr (ip) z JO 

In case (A), AA' 1U ^ 0, and both, Eq. (|A2|) and (|A3|> . are contributing to Tii.b, leading to a cancellation of 
the contrubution from the change in the projected meniscus area of Eq. I|A3(1 by the second term in Eq. I|A3(1 . 
Inserting the decomposition of f{p>) from Eq. (J3J) (omitting the colloid label i) we find 

n h ^ J /"d? [/(</>) - /i] 2 



2 

7T7 

T 







2(P -M 2 + 4 V |P„ 



m>l 



(A4) 



For the case (B3) there are also tilt fluctuations of the vertical axis of the three phase contact line circle which 
can be parametrized by a boundary height according to f — P\e lip + P-ie~ lip , the projection of the tilted circle 
onto the reference plane z = is an ellipse, and we find for the boundary Hamiltonian according to Eq. l|A3jl 

n h ~ ^ 7 (|Pi| 2 + |P-i| 2 ) (as) 

Note, however, that we can write 7ib,i as in Eq. Q for all cases for the boundary conditions, (A) and (B1)-(B3), 
since the integration measure Vfi for the contact line is constructed such that only the actual terms for the 
respective cases contribute, see Sec. Ill Al 

APPENDIX B: EXPANSION OF THE GREENS FUNCTION 

In this appendix we derive the multipole expansion of the Greens function G(|x|) w — (l/27r) ln(7 e |x|/2A c ) 
between two "charged" (charges generating the auxiliary field ifti) and isolated regions (dSi trc { in Sec. IIII Al and 
Si.ref in Sec. IIIID|) . The method we apply is known from electrostatics and goes back to Schwinger (c.f. 44]) 
and is referred to as the Schwinger technique in the literature. In doing so, we use the analogy of our problem to 
two-dimensional electrostatics. As starting point we use the fact that the logarithm is the generating function 
of the Gegenbauer polynomials. Using r — r'| = \/r 2 + r' 2 — 2rr' cos(p — p') and exploiting some properties of 
the Gegenbauer polynomials [2^ we find 

In |r -r'| = lnr + ^ y (—} cos(l<p - lip') , (Bl) 



I 

l>i 



where we assumed r' = |r'| < r = |r|. Comparing Eq. 1|B1|) to the Taylor expansion of the logarithm, we find 

ti lnr ~ | (^y ((Mto-w +e -i{i<p-iv')^ i > o • (B2) 

On the other hand, introducing £± = d x ± id y and using lnr = lnr = 0, we can write 

(-*•'■ V)< ,__ (-1)' ' ' 



II I! 
Identifying Eqs. I|B1J1 and l|B3|l we obtain 



lnr=^— J- - g. +e u * £. lnr (B3) 



£ ± lnr = —(—2) (I — 1)!— j — (B4) 
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With Eq. (|B4(1 at hand and for xi = ri and X2 = d + r2 residing on different circles we obtain for the Greens 
function G(|x 2 — Xi|) 



J_ , / 7 c |d + r 2 - ri 
2tt n V 2A r 



2?r V 2A r 



2?r V 2A r 



1 E 

2tt ^ 



-ri • V)' 1 (r 2 • V)' 2 



2vr ^ /i! 

11^2=0 
i!+i 2 >l 



In r 



2tt .4^. ^i!Z 2 !2' 1+ ' 2 



+ 



lnr 



z 1 ,z 2 =o 

'l+'2>l 



2tt V 2A C 



2^ ^ Zx + ^V ii ^d'^ 



g-i(/iyi+/2V2) _|_ g»(/l.V>l+i3V>2) 



(B5) 



In Sec. 1111 Al we also need the multipole expansion of G(|x 2 — Xi|) for xj residing on the circumference dSi <Ie { of 
the same circle in order to calculate the elements of the self-energy matrix, Eq. i|17|) . Using Eq. (|B1[1 , we obtain 
with |xi — xa| = r \/2 - 2cos(<^i - </? 2 ) 



G(| Xl -x 2 |) ~ -—In 
= -i-ln 



2tt V 2A. 



7c 
2A C 

J c ro 



Z7T 



2vr ^ 2? 

2>i 



(B6) 



Here the prerequisite r' < r of Eq. I|B1(I is not fulfilled, and the series in Eq. (|B6J) is not convergent. It has to be 
understood in a formal sense, since it only provides the Fourier coefficients for a finite number of modes which 
actually contribute to the effective interaction in the long-range regime (see Sec. MI Al for details). 



APPENDIX C: KARDAR'S METHOD: CALCULATION OF THE SELF-ENERGY 



The self-energy part, Eq. (|49|l . is evaluated quite similarly as in Ref. [T^. We eliminate the 8 functions by 
introducing conjugate multipole moments of the auxiliary fields: 



d 2 x(r/r ) lml e- lm *TP{x) = / # im exp i$ 



S I * im - 

This brings iJi, sc if = cxp{— ksT/(2j) H.i. sc if } into the form 

Zi.seif = J Y[d^ m J Vi>i exp | -M^ J d 2 x J 



d 2 x(r/r ) |m| e- lmv V(x) 



(CI) 



, . d 2 x / d 2 x'^(x)G(|x-x'|)^(x) 

27 Js lirof Js,,^ 



/ d 2 x^(x)V - (P im + * im )e im ^ + i V (P im + ^ m )tf 



(C2) 



The functional integral / P^j in Eq. I|C2|) can be converted into a functional integral over a constrained height 
field h(x); this corresponds to a reversion of the step from Eq. I|14fl to Eq. Q45[l. 



-^i.self = / ]~J<f* im exp{« ^ (Pirn + *im)*im} / ^ J| <5 I 7l(x;) - ^ 



x exp 



2fc B T 



7 'd 2 * 



(v/,) 2 + 



(C3) 
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In Eq. (|C3|) . the 5- functions describe the pinning of the field h in the region Si tTe {. This contribution can be 
evaluated directly, such that the remaining functional integral reads 



-2;, self 



A c — >oc 



exp 



2-7T7 



m>l 



x exp ■ 



d 2 x 



A 2. 



(C4) 



where the (5-functions fix /i(x) at the boundaries dSi jTe { of the integration domain. Splitting the auxiliary 
field into two parts as in Eq. ©, h = h a + hi , where (—A + A f T 2 )/i = with the boundary conditions 

ho{xi)\dS i>rt>f = Em=-A + *im) e! ' nw and ^i(xi)|aSi, r ef = 0, and applying Gauss' theorem to the integral 
in the exponent of Eq. (|C4|) leads to 



Zi,seU = / II *™ ex P S ■ 



27T7 



X] IHI-Pim + %m\ 2 + i (Pirn + *im)% 



m>l 



x exp 



<p dxh (x.)\7h (x) 



2k B T, 

fvhi J[ 5(/n(x))exp. 

« ^- d o 



2 7 7[R2\s iircf 



d 2 s 



A? 



(C5) 



The functional integral over h\ only yields a constant factor independent of any multipole moment, which will 
be disregarded below. To compute the line integral in Eq. I|C5|1 . we write the general solution of the differential 
equation for ho in R 2 \Si, r ef as (c.f. Sec. IIII Btl /io(x) = Y^nS-Kmi? I '^c)/K m (r / \ c ))C m e lmip . By comparison to 
the boundary conditions the coefficients are determined straightforwardly as C m = Pim + ^im- Then, the line 
integral evaluates to 2^\P im + ^i m \ 2 f(m) with f(m) = \m\ (\m\ > 1) and /(0) = — l/(ln(7 e ro/ 2A C )), such that 
the self energy part reads 



^i,sclf 



, exp 



47T7 
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and thus Tii^oU is given by (up to unimportant additive constants) 

■ 1 

2tt/(to) 



^2, self / ^ 



(C6) 



(C7) 



m>l 



[1] J. D. Joannopoulos, Nature 414, 257 (2001). 
[2] P. Pieranski, Phys. Rev. Lett. 45, 569 (1980). 

[3] F. Ghezzi and J. C. Earnshaw, J. Phys: Condens. Matt. 9, L517 (1997). 

[4] F. Ghezzi, J. C. Earnshaw, M. Finnis, and M. McCluney, J. Colloid Interface Sci. 238, 433 (2001). 
[5] J. Ruiz-Garcia and B. I. Ivlev, Mol. Phys., 95, 371 (1998). 

[6] J. Ruiz-Garcia, R. Gamez-Corrales, and B. I. Ivlev, Phys. Rev. E 58, 660 (1998). 
[7] D. Stamou, C. Duschl, and D. Johannsmann, Phys. Rev. E 62, 5263 (2000). 

[8] M. Quesada-Perez, A. Moncho-Jorda, F. Martinez-Lopez, and R. Hidalgo- Alvarez, J. Chem. Phys. 115, 10897 
(2001). 

[9] R. P. Sear, S.-W. Chung, G. Markovich, W. M. Gelbart, and J. R. Heath, Phys. Rev. E 59, R6255 (1999). 
[10] W. Chen, S. Tan, T.-K. Ng, W. T. Ford, and P. Tong, Phys. Rev. Lett. 95, 218301 (2005). 

[11] M. G. Nikolaides, A. R. Bausch, M. F. Hsu, A. D. Dinsmore, M. P. Brenner, C. Gay, and D. A. Weitz, Nature 424, 
1014 (2003). 



24 



[12] A. Wiirger, and L. Foret, J. Phys. Chem. B 109, 16435 (2005). 

[13] M. Oettel, A. Dominguez and S. Dietrich, Phys. Rev. E 71, 051401 (2005). 

[14] M. Oettel, A. Domiguez and S. Dietrich, J. Phys.: Condens. Matter 17, L337 (2005). 

[15] M. Oettel, J. Phys: Condens. Matter 17, 429 (2005). 

[16] R. A. McGough, and M. D. Miller, Phys. Rev. A 34, 457 (1986). 

[17] A. Kovalenko, and F. Hirata, Phys. Chem. Phys. Chem. 7, 1785 (2005). 

[18] D. Jasnov, in Phase Transitions and Critical Phenomena 10, ed. by C. Domb and J. L. Lebowitz (Academic, London, 
1986). 

[19] R. Golestanian, M. Goulian and M. Kardar Phys. Rev. E54, 6725 (1996). 

[20] M. Goulian, R. Bruinsma, and P. Pincus, Europhys. Lett. 22, 145 (1993); 23, 155(E) (1993). 

[21] H. Lehle, M. Oettel, and S. Dietrich, Europhys. Lett. 75, 174 (2006). 

[22] Bordag M., Mohideen U. and Mostepanenko V. M., Phys. Rep. 353(2001). 

[23] M. E. Fisher and P.-G. de Gennes, C. R. Acad. Sci. Paris Ser. B 287, 207 (1978). 

[24] M. Krech, The Casimir Effect in Critical Systems , (World Scientific, Singapore, 1994) 

[25] A. Ajdari, L. Peliti, and J. Prost, Phys. Rev. Lett. 66, 1481 (1991). 

[26] H. Li and M. Kardar, Phys. Rev. A 46, 6490 (1992). 

[27] Statistical Mechanics of Membranes and Surfaces, Proceedings of the Fifth Jerusalem Winter School, edited by 

D. R. Nelson, T. Piran, and S. Weinberg (World Scientific, Singapore, 1989). 
[28] S. K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997). 

[29] T. Emig, R. L. Jaffe, M. Kardar, and A. Scardicchio, Phys. Rev. Lett. 96, 080403 (2006). 

[30] R. L. Jaffe and A. Scardicchio, Phys. Rev. Lett., 92, 070402 (2004). 

[31] B. Derjaguin, Kolloid-Zeitschrift 69, 155 (1934). 

[32] H. Gies, and K. Klingmuller, Phys. Rev. Lett. 96, 220401 (2006). 

[33] R. Buscher and T. Emig, Phys. Rev. A 69, 062101 (2004). 

[34] R. Golestanian, Phys. Rev. E62, 5242 (2000). 

[35] S. Safran, Statistical Thermodynamics of Surfaces, Interfaces, and Membranes (Westview Press, Boulder, 1994). 
[36] H. Li and M. Kardar, Phys. Rev. Lett. 67, 3275 (1991). 

[37] H. Kleinert, Path Integrals in Quantum Mechanics, Statistics and Polymer Physics (2nd edn.), (World Scientific, 
Singapore, 1995) 

[38] P. A. Kralchevsky, V. N. Paunov, I. B. Ivanov, and K.Nagayama, J. C oll. Int. Sc. 151 (1992), 79. 
[39] C. Bachas, P. Le Doussal, and K. J. Wiese, preprint (hep-th/0606247). 
[40] H. Kaidi, T. Bickel and K. Benhamou, Europhys. Lett. 69, 15 (2005). 

[41] W. H. Press, S. A. Teukolsky, W. T. Vctterling, and B. P. Flannery, Numerical Recipes in C++ (2nd edn.), (Cam- 
bridge University Press, Cambridge, 2002). 
[42] L. Helden, private communication. 

[43] D. G. A. 1. Aarts, M. Schmidt, and H. N. W. Lekkerkerker, Science 304, 847 (2004). 

[44] J. Schwinger, L. L. DeRaad, Jr. Kimball A. Milton, Wu-yang Tsai, Classical Electrodynamics, Perseus Books, 
Reading (1998). 

[45] M. Abramowitz and I. A. Stegun (Eds.), Handbook of Mathematical Functions (Dover, New York, 1974). 
[46] Since Pi m = P*_ m , real and imaginary parts of Pi± m are not independent of each other, and therefor we define the 
measure dP irn dPi- m = dReP irn dImPi_ m . 



